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Abstract 

We present a quantitative analysis of the electromagnetic pion form factor in the 
light-cone sum rule approach, including radiative corrections and higher-twist effects. 
The comparison to the existing data favors the asymptotic profile of the pion distribu- 
tion amplitude and allows to estimate the deviation: (J du/u(j)Tj{u))/{f du/ucp^^u)) = 
l.lib 0.1 at the scale 1 GeV. Special attention is payed to the precise definition and 
interplay of soft and hard contributions at intermediate momentum transfer, and 
to matching of the sum rule to the perturbative QCD prediction. We observe a 
strong numerical cancellation between the soft (end point) contribution and power 
suppressed hard contributions of higher twist, so that the total nonperturbative cor- 
rection to the usual pQCD result turns out to be of order 30% for ~ 1 GeV^. 
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1 Introduction 



There is a clear tendency for QCD-oriented experimental studies to go for more and more 
exclusive channels. All future plans also call for very high luminosity and would therefore be 
perfectly suited for the investigation of exclusive and semi-exclusive reactions. A problem 
which hinders all attempts to implement these projects is the lack of truly quantitative 
QCD predictions. It is widely anticipated, see e.g. [|^, 0, H, that for experimentally 
accessible values of the momentum transfer, the perturbative QCD factorization for hard 
exclusive reactions receives non- negligible corrections from the so-called soft, or end- 
point contributions, which are essentially nonperturbative. One practical difficulty is that 
soft corrections can in many cases be mimicked (numerically) by modifying the shape of 
hadron distribution amplitudes. An agreement of perturbative predictions with the data 
cannot, therefore, be used to claim smallness of end-point effects which have to be estimated 
independently using a certain nonperturbative approach. Creating a systematic framework 
for a study of soft end-point corrections is becoming, thus, increasingly timely. 

It has been suggested that the soft end-point contribution to the pion electromagnetic 
form factor can be estimated in a largely model-independent way within the framework 
of light-cone sum rules [0]. The aim of the present paper is to put this technique on 
a more quantitative footing. To this end we calculate the radiative correction to the 
light-cone sum rule, elaborate on the scale dependence and demonstrate how the sum 
rule estimates of the end-point effects can naturally be combined with the NLO QCD 
perturbative calculation. In addition, we estimate the twist 6 contribution to the sum rule 
due to the quark condensate and find this correction to be small. 

The presentation is organized as follows. In Sect. 2 we remind basic ideas of the 
light-cone sum rule approach and derive the simplest sum rule. Sect. 3 is devoted to the 
calculation of the radiative correction and separation of soft and hard effects. As expected, 
we find that in the — >■ oo limit the form factor is dominated by the hard rescattering 
contribution alone, while to the accuracy both soft and hard contributions have to 

be taken into account. Higher-twist corrections to the light-cone sum rule are considered 
in Sect. 4, while Sect. 5 contains the results of our numerical analysis. Matching of the 
sum rule with the NLO perturbative predictions is discussed in Sect. 6. Finally, in Sect. 7 
we summarize. The paper contains two appendices where we collect some useful but bulky 
expressions and present the relevant formulae for light-cone distributions of the pion. 



2 The method of hght-cone sum rules 

The approach is based on the study of the correlation function 

where j'^ = d'j^'y^u and j^™ = euU'~^uU + edd'y^d is the quark electromagnetic current. With 
= and = —q^ fixed, the correlation function (|l|) depends on a single invariant 
variable s = {p — qY- The contribution of the pion intermediate state equals 

T^,(p, q) = 2iU{p - q)^p^F^[Q^)— } , (2) 
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Figure 1: The tree- level contribution to the correlation function in Eq. (|l]). 

where fjr is the pion decay constant and F-j^iQ"^) is the pion electromagnetic form factor. On 
the other hand, at large negative {p — qY ^"^^^ the correlation function can be calculated 
in QCD, in full analogy with the 7*7*7? transition form factor. A common idea of all QCD 
sum rules is a matching between the QCD calculation at Euclidean momenta and the 
dispersion relation in terms of contributions of hadronic states, which allows to estimate 
the hadronic quantity of interest. Specifics of the light-cone sum rules is how exactly the 
QCD calculation and matching are done. To illustrate this point, consider the contribution 
of the simplest diagram in Fig. 1: 



T, 



27r2 



d X— (0|[e„(i(0)7^ ^'-iy'-i^u{x) - edd{x)-f,y ^ltil5u{0)]\'iT^{p)) . (3) 



X 



Expansion of the remaining nonlocal matrix elements around the middle point in a formal 
Taylor series generates the Wilson operator-product expansion in contributions of local 
operators of increasing dimension 



O 



tl,fJ.l,fl2,...,tJ.n 



(4) 



Df^i = — igA^ being the covariant derivative and = — D^- Restricting ourselves 
for the moment to operators of the lowest twist (highest Lorentz spin) we consider the 
relevant reduced matrix elements 



X 



Ml 



• • •^.„(0|O;,,^,,,,...,,j7r+(p)) = ^/.p,(px)"((0„)) + . . . . 



(5) 



They are related, as first found in H, to the moments of the pion distribution amplitude 



(0|ci(0)7^75u(x)|7r+(p)) = ip^U j dMe-*"^>^(M,/x^ ~ x"^) 





{{On)) 



du (1 - 2uY^M 



(6) 
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Substituting Eq. in the expansion of and integrating over we obtain for the 
contribution of Fig. l[| 

+ other Lorentz structures, (7) 
where s = {p — q)"^. To construct a sum rule, we make the Borel transformation 

> exp[-ml/M^] , 



1 



ml — s 



1 



exp[-gV^^ 



(Q2 _ s)n (M2)"-ir(n) 

introducing a new variable (the Borel parameter), and equating the Borel-transformed 
versions of Eqs. (H) and (^. For simplicity we neglect the continuum subtraction here. 
Neglecting the pion mass, the result reads 

n-l\ 1 f-2Q 



n=2,4,... fc=l 



This sum rule is, however, completely unsatisfactory! 

Indeed, QCD sum rules are generally expected to hold in a certain interval of values of 
the Borel parameter, such that contributions of both higher resonances and higher orders 
of the OPE are simultaneously suppressed. It is easy to see that in the present situation 
these two conditions are contradictory, unless is sufficiently small. Indeed, on the one 
hand, one has to keep small, of order 1 — 2 GeV^, to suppress the contribution of, e.g., 
the Oi-meson intermediate state. On the other hand, for a fixed the higher order terms 
on the r.h.s. of the sum rule are enhanced by factors [Q"^)^ and for > the OPE 
expansion breaks down. 

An escape suggested in is to avoid the Wilson short-distance expansion altogether 
and write the answer for the diagram in Fig. 1 directly in terms of the pion distribution 
amplitude. The expansion parameter then becomes the twist of the operators rather than 
their dimension. Using Eq. (^ and the definition of the pion distribution amplitude in 
Eq. (ID we obtain to leading twist accuracy, instead of Eq. (^ , a compact expression 

1 

T^u = 2tUp,p, [ du^^^ + ..., (10) 
J uQ"^ — us 



where m = 1 — u . Making, once again, the Borel transformation, we get the simplest 
light- cone sum rule ^ 

FAQ') = J duifM^W (-^) • (11) 



^The terms with odd n vanish because of G-parity. 
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This sum rule is perfectly well behaved at — oo and it is instructive to trace how 
the above-mentioned difficulties of the standard approach have been resolved. Because 
of the strong exponential suppression factor, the important region of integration over the 
momentum fraction variable u gets shifted, in the large-Q^ limit, to the end-point region 
1 — M ~ M^/Q"^. In this regime, the virtuality of the quark (the denominator in Eq. ([To|) ) 
remains all the time of order M^, as —>■ oo. The deficiency of the short-distance 
expansion is now clearly seen as originating from the wrong expansion parameter {Q'^ — s)/2 
(cf. Eq. (|^)) corresponding, effectively, to the expansion around the symmetric point 
« = 1/2 

To be somewhat more quantitative, we have to make the usual continuum subtraction. 
This is trivial in the case at hand, since the expression ([T0|) is easily converted to the form 
of a dispersion integral over s = {p — q^. All we have to do is to truncate this integral at 
a certain threshold Sq, called the interval of duality. The result ||^ is that the integration 
over the momentum fraction is cut from below at the value 

wo = QVK + Q'). (12) 

In addition, the pion distribution amplitude has to be taken at the scale corresponding to 
the quark virtuality 

fil = uQ^ + uM^ . (13) 

Implementing these small improvements, we obtain the leading-twist leading-order light- 
cone sum rule @] 

1 



(14) 



wo 



The crucial advantage of the light-cone sum rule approach is that it allows to incorporate 
the information on the end-point behavior of the pion distribution amplitude ipn{u) 
1 — u. In the limit ^ oo the integration region in Eq. ( [T^ ) shrinks to a point m = 1 so 
that one obtains 

J, (15) 



where (p'^{0) = {d/du){pT,{u)\u^o = ^V^^l)- The Borel variable corresponds to the 
(inverse) distance at which the matching is done between the parton and hadron represen- 
tations. 

The expressions in Eqs. ([T^), (|TB|) present a typical 'soft' or 'end-point' contribution to 
the pion form factor which is sensitive to the pion wave function at a low normalization 
point and comes from large transverse distances of order b ~ Sq^^"^- 

To illustrate this point, write the four- dimensional integration in Eq. (|I]) as a product of 
two two-dimensional integrations in longitudinal and transverse (to p and q) coordinates. 



similar deficiency of tlie sliort-distance expansion in the case of heavy-to-liglit correlation functions 
is demonstrated in M. 
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Figure 2: The transverse-distance separation between the quark and the antiquark in the leading- 
order hght-cone sum rule ( 14 ) in the large limit for typical values of the sum rule parameters 
so = 0.7 GeV^ and = 1.0 GeV^. 



Leaving the transverse integration intact, a short calculation gives for the r.h.s. of Eq. (|TT 



d b J du ifniu 





I 

47r 



exp 



4 mM2 J 



(16) 



The distribution of transverse distances in the diagram in Fig. 1 is, thus, gaussian, with 
the average transverse size (6^) = 4/{uM'^) controlled by the value of the Borel parameter. 
One also sees that the scale of the distribution amplitude in Eqs. (|13|), (p^ is determined 
by the weighted average of the momentum transfer and the (inverse) transverse distance 
between the quarks, as expected on general grounds [ p!0|] . 

Including the continuum subtraction modifies this distribution rather significantly as 
the small-6 region is dominated by high-mass excitations and gets suppressed. After some 
algebra we obtain the sum rule equivalent to Eq. ( p!^ but with an explicit separation of 
different transverse distances: 



1 

Air 



d% / du ip-n-iu) e 



USQ—UQ 

-uQyiuA4^) f ^-t/{uM^) 



UO 







MVbH) 



47rQ4 



so 



d b dse 



dt{s-t)Jo{VbH) 



(17) 



where Jq is the Bessel function. The resulting transverse-distance distribution (normalized 
to unity at 6 = 0) is shown in Fig. 2. The dependence on both and the Borel parameter 
is actually very weak and the overall scale of transverse distances is determined almost 
entirely by the value of the continuum threshold. Because of this, for ^ sq the pion 
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Figure 3: The radiative corrections to the correlation function in Eq. (|^). Dashed hues denote 
virtual gluons. 



distribution amplitude has to be taken at the scale /x^ ~ Sq, rather than at yU^ ~ f\. 
The width of the 6^-distribution in Fig. 2 should be compared with the electromagnetic 
pion diameter squared (2_R^™)^ ~ 2 fm^. 



3 Radiative corrections 



3.1 General case 



In order to improve the accuracy of the light-cone sum rule (p!^, one has to calculate 
the 0{as) radiative corrections to the leading-order correlation function (0). The cor- 
responding Feynman diagrams are shown in Fig. 3. The calculation is straightforward, 
albeit tedious, and technically similar to the calculation of the radiative correction to the 



7*7*7r transition form factor for different photon virtualities |]T2[. We handle ultraviolet 
and infrared coUinear divergences by dimensional regularization in the MS scheme. Due 
to the fact that the diagrams contain two 75 matrices (one from the axial-vector vertex 
and one from the pion projection) there is no 75 ambiguity. We have also checked that the 
collinear divergences are absorbed in the definition of the scale-dependent pion distribution 
amplitude ip-,,- Our result for the twist 2 part of the correlation function (|I]) to 0(0;^) can 
be represented in a form of convolution of (f-^ with the hard scattering amplitude 



1 





■^It can be shown that this change of scale takes into account the continuum subtraction in the running 
couphng, cf |jn|. 



where the leading-order result was already given in Eq. (|TI1|): 



u 



uQ'^ — us ' 

and the radiative correction to the hard scattering amplitude equals 



uu{Q'^ + sYiuQ'^ — us) 



- 9uu{Q^ + sY 



(19) 



- \Q%(3u - 2) + Q2s(5 _ Quu) + s\(3u - 2)1 In , 



uQ — us 



+ [Q'^uu + Q^s{l + 2uu) + s'^uu] In^ 



uQ"^ — us 



+ u \-2Q% + 5Q'^s + 5^(1 + 2u)] In — - us [Q^fl + u) + su] In^ — 



+ u [Q\l + 2u) + SQ^s - 25^^] In - mQ^ [g2^ ^ ^ ^^^j jj^2 v 



.(20) 



To obtain the radiative correction to the light-cone sum rule for F-^, one has to calculate the 
imaginary part of Hi in the variable s. The resulting expression is presented in Appendix 
A. After continuum subtraction and Borel transformation, the sum rule reads: 





where 

J^l(«,M^.o) = 

-|2 



Q{u - no) J-S(w, So) + Q{uo - u)J'2diu, M\ s^) 



, (21) 



exp 



mil 



uM^ / r 47r ^ 



-9 



TT 



3 In ^ + 3 In — ^ - In^ ^ - In^ 



u^^ 



uQ^ 

UIJ? 



so 



(is Q e 



u{Q^ + s) 



5s + 1 + 2 In ^ + 2 ^ + s In 



u 



uQ^/u 

2Q^fQ^ + s 2M2 + Q2 + 5 _p\ _p 



+ 



M \ S 



+ 



(is Q'^e 



2„-s/M2 







uu{Q^ + s) 



2u{Q^ ~s + s\n^] + ( -g2 + 5s + 2(g2_s)in4 



^^^' + ^)^ -3 + 21n4nin^ 



M2 



+ 2^e-«/^^^ln^ln 



Mo 



(22) 



and 



47r" 



so 



ds Q'^e 



-s/M'^ 



uiQ^ + s) 



2 - s + s In ■ 



+ 



u 



g2 + 5s + 2(Q2_s) In. 



s{Q^ + s) 
M2 



21n4 



iu4 



UU 



21ui| 



(23) 



Here p = uQ"^ — us and po = uQ'^ — usq = (1 — u/uq)Q'^. The superscript '*^2)' indicates 
the leading twist 2 contribution. Higher twist terms will be added in the next section. 



We interpret the parts of Eq. with ^hard('"' ■^o) and M^, sq) as "hard" and 

"soft" contributions to the pion form factor, respectively, defined with the explicit cutoff 
in the momentum fraction m = mq ~ 1 — Sq/Q^. This separation will be discussed in detail 
below. 



:(2) 



3.2 Study case: asymptotic distribution amplitude 



For the asymptotic shape of the pion distribution amplitude ^^{u) = 6u{l — u) the 
momentum-fraction integration in the radiative correction can easily be done analytically, 
with the simple result 



so 



[s + Q 



2\4: 



1 + 



a.C. 



7T 

y 



6-ln^^ + A + ^ 



(24) 



All the scale-dependent logarithmic terms cancel in this expected. For large 

^ So one can expand the sum rule (|2^) in powers of l/Q^- 




so 



dsse~'/''Ul 



Arc 



3 s 



0(l/g^). (25) 



so 



To interpret the leading term, we notice that the integral J ds e~^^^ can be related to the 



pion decay constant through the QCD sum rule |]13| 



12M2 



+ 



81M4 



Txas{qqY + 



(26) 
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The perturbative correction and the gluon- and quark-condensate contributions involve 
an extra power of and are absent, therefore, in our approximation. Substituting 



/ ds e" 





47r we obtain 



(27) 



which coincides with the classical result [^]. 

It is easy to see that the 0{1/Q^) contribution to the form factor comes entirely from 
the term which we have identified as "hard" , while all power suppressed corrections involve 
both hard and soft contributions. In particular, we obtain to 0{1/Q'^) accuracy 



^ as, hard 

TT 



Uo 



du if 



So) 



2ttQ^ 



so 



dse-^/''Ul- 



2^ + ln4+f3-21n4 
s fi^ \ /i^ 



In- 



So 



(28) 



The soft part is then just what is left when this hard contribution is subtracted from the 
total result in Eq. (pSl) . Notice that the separation of soft and hard contributions depends 
on the collinear factorization scale, even for the asymptotic distribution amplitude. We 
will elaborate on this dependence in what follows. 
In the local duality limit —>■ oo one obtains 



71 as, hard 



27rQ2 



Soil 



So_ 



and 



13 
Y 



So So 



So 



So 



(29) 



^as,soft|-Q2'j _ 3Sq ^ SasCp ^2 



- + In^ — 

2 So 



ln2^ + 21n^ + 31n^ 

fJ- So So 



(30) 



where Sq ~ 47r^/^, cf. Eq. (pGf). Note that the 0(1/Q^) hard contribution is large and 
negative, while the soft radiative correction O^ag/Q"^) is positive, unless ^ So,/i^. This 
implies considerable cancellations in the sum of the soft and hard contributions so that in 
order to make their separation physically meaningful one must assume a low value of the 
factorization scale /i^ ~ so.Q 

Finally, notice the double-logarithmic contribution ~ \n^Q^/s in Eq. (|2^) which is 
reminiscent of the Sudakov logarithms discussed in []T0| . A typical size of these corrections 
is of order In^ Q'^/ sq which for sq ~ 0.7 — 0.8 GeV^ and ~ 1 — 10 GeV^ is much less than 
In^ Q^/Aqqj-, with Aqcd ~ 200 MeV, as usually assumed. For this reason, exponentiation 
of Sudakov corrections is numerically not important in the present approach. 

"'it is easy to see that for /z^ = there are double-logarithmic contributions ^ ln^(Q^/so) to Eqs. ( [29| ) 
and (pO) which have opposite sign and partially cancel in the sum. 
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3.3 The expansion in general case 

For a generic pion distribution amplitude the NLO light-cone sum rule (0) again simplifies 
considerably upon the expansion in powers of l/Q^. We obtain to 0(1/ Q^) accuracy 



a 



271 



du 



u 



so 



9 + ^vr^ + In 4^ - In^ 
3 iJ,^ 



s 



as fdsse ''/^^^ 





(2 In — - 3) / du 



iPjr{u) -U(p'{0) 



+ (2 In — - 8) Idu 



u 



(31) 



The scale dependence cancels to the required accuracy, since [|4 
d 



a. 



-y,;(0,/i) = ^ 

dm 11 71 



du 



u 



(32) 



The contribution of hard rescattering equals 



27r' 



*0 ,,,9 1 

dse-''^' 



du 



(p^{u) 



u 



as ^ fds se '^/^^ 





(2 In — - 3) / du 



+ (21n. 



du 



U 



so 



dse-^ 



s ( 2 In — - 3 I In 



So - s 



2so } , (33) 



and the soft contribution is identified as the difference F^°^*(Q^) = Fj,{Q'^) — F^'^'^'^{Q'^). 
Note the term proportional to V'^(O) = —ip'^{l) in the last line of Eq. (^) which is concen- 
trated at the end-point but enters as part of the hard contribution. 
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A few comments are in order concerning this expansion. First, consider the leading 
asymptotic 0(1/Q^) term. Substituting, as above, j^^ds e~*/^^ — >• 47r^/^, this contribution 
can be rewritten as 



(34) 



where we used that j^dvip^iy) /v = 3. This expression is similar, but does not yet coincide 
with the perturbative QCD result H] 



2^ _ 87raJ,^ 



1 2 

du"^ 



u 



(35) 



It is easy to convince oneself that the missing corrections to the asymptotic pion distribution 
amplitude in the first integral in Eq. (|3^ as well as the missing nonperturbative corrections 
to are supplied by eventual higher-order and higher-twist corrections to the sum rule. 
Since such corrections are difficult to evaluate directly, one may try to improve the light- 
cone sum rule by combining it with the known full NLO perturbative calculation. Such a 
possibility will be discussed in Section 6. 

Second, the structure of the 0(1/(5^) power correction to the pion form factor is very 
similar to the heavy quark limit of the light-cone sum rule for B — >■ 'kcv decay considered 
in In particular, note the l/v? weight factor in the integral over the pion distribution 
amplitude, the structure of double-logarithms and, finally, the cancellation of the coUinear 
factorization scale-dependence by the same mechanism. In both cases, the distinction 
between soft and hard contributions necessitates a kind of generalized "plus-distribution" 
subtraction of divergent integrals over the pion distribution amplitude at u — > 1, as it 
is done in the second line of Eq. (^31). This feature seems to be general, whereas the 
distribution of finite terms ~ </'^(0) between the hard and the soft contributions is arbitrary. 
The expression for such terms in the last line in Eq. (|33D corresponds to the particular 
definition ( PT| ) with a rigid cutoff in momentum fraction. Although such a definition is the 
most intuitive one, it is not unique and, as seen from Eqs. (pi]), (p3D, introduces rather 
cumbersome "surface terms" ~ V^^(O) which appear both in hard and soft contributions 



and cancel in their sum. An interesting alternative ||T^ which we do not pursue in detail 
in this work is to define the separation between hard and soft contributions order by order 
in the expansion using "plus-distributions" to regularize the divergent momentum- 

fraction integrals. To the accuracy of Eq. (pT]), this procedure corresponds to the definition 
of the hard contribution as given in the first three lines in Eq. (|HB| ) omitting the "surface 
term" . The soft contribution is given then by the second line in Eq. (|3l| ) . 

Last but not least, having in mind that the separation of hard and soft contributions is 
ambiguous, one may add them together and consider their sum as a 'total nonperturbative' 
power correction to each order in the expansion. Inspection of Eqs. ([25|), (P^), ( pUD 
suggests that soft and hard corrections in general have opposite signs and partially cancel 
in the sum. We postpone the detailed discussion of this issue to Sect. 6 where we summarize 
our numerical results. 
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Figure 4: The relative contributions of hard (dashed) and soft (dotted) regions in the leading- 
twist light-cone sum rule ( ^T|) for the pion form factor (the sum: solid lines). The upper figures 
correspond to asymptotic and the lower figures to CZ pion distribution amplitude. The three 
figures from left to right correspond to the choice of factorization scale = Q^-, fJ^ = and 
= uQ"^ + uM'^, respectively. 

3.4 Numerical estimates 

Results of the numerical evaluation of the sum rule (pT]) are shown in Fig. 4 by solid curves, 
for So = 0.7 GeV^ and for a typical value of the Borel parameter = 1.0 GeV^. (The 
choice of input parameters is discussed below in Sect. 5). Soft and hard contributions 
are shown by dotted and dashed curves, respectively. The results are plotted using the 
asymptotic ^p'^{u) = Qu{l — u) and the Chernyak-Zhitnitsky (CZ) ip'^'^{u,fi = 1 GeV) = 
30 MM (2m — 1)^ pion distribution amplitudes, for three different choices of the factorization 
scale: /i^ = Q^, /i^ = sq (see the discussion above) and /i^ = /z^ = uQ^ + uM^, according 
to Eq. (|13D 0. In all calculations in this paper we use the two-loop QCD running coupling 
with = 336 MeV corresponding to Osil GeV) = 0.48 and agiso) = 0.59. It is seen that 
hard contribution to the form factor defined with a "natural" momentum-fraction cutoff 

'"^ If the momentum-fraction dependent scale Hu is used, it is implied that as{i-iu) is inserted inside the 
M-integrals. 
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Figure 5: Average momentum fraction m as a function of for the asymptotic (solid) and CZ 
(dashed) distribution amphtudes at = 1.0 GeV^. The scale is /i^ = (1 — u)Q'^ + uM"^ . 

remains small and negative for the main part of the interesting region of Q^. 

Furthermore, in Fig. 5 we show the average value of the momentum fraction u in the 
integral in Eq. (|2l|) calculated as a function of for the asymptotic (solid curve) and 
CZ (dashed curve) distribution amplitudes. This average value turns out to be very large, 
and, contrary to usual expectations, does not depend significantly on the shape of the pion 
distribution amplitude. 

Negative contribution of the hard-rescattering mechanism may appear unexpected and 
counterintuitive. We emphasize, however, that the separation between hard and soft terms 
is ambiguous and depends on their definition — this is, in fact, the main lesson to be learnt 
from our analysis. Note that the scale dependence is much more pronounced for hard and 
soft contributions taken separately than for their sum. 



4 Higher-twist corrections 
4.1 Twist 4 contributions 

The operator-product expansion of the correlation function (|I]) near the light-cone a;^ = 
can be continued beyond the leading twist 2 approximation ([T0|). This procedure yields 
higher- twist corrections to the light-cone sum rule (^TJ). They are suppressed by additional 
inverse powers of and . Physically, the higher-order terms of the light-cone expansion 
take into account both the transverse momentum of the quark-antiquark state and the 



contributions of higher Fock states in the pion wave function. As explained in [^], these 
two effects are indistinguishable due to QCD equations of motion . 

Next to the leading twist 2 term, the correlation function (|1|) receives several twist 4 
contributions ^ First of all, one has to take into account the twist 4 components of the 
quark-antiquark matrix element (0|(i(0)7^75M(a;)|7r"'"(p)) in the diagram in Fig. 1. Further- 

^Twist 3 contributions to dl^) are proportional to and vanish in the chiral hmit adopted here. 
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a) b) c) 

Figure 6: The light-cone expansion of the quark propagator in the correlation function 



more, the gluon emission from the virtual quark should be included yielding the diagram 
of Fig. 6a with the twist 4 quark- ant iquark-gluon distribution amplitudes of the pion. To 
calculate this diagram, one makes use of the light-cone expansion of the quark propagator 
p!6| given in Appendix A. The definitions of all relevant twist 4 two- and three-particle 



distribution amplitudes [|T5|, are collected in Appendix B. 

The corresponding calculation has been carried out in 0]. The twist 4 contribution to 
the correlation function can be written in the following compact form 



TifJ = 2zp^pJ^ I , (36) 



{uQ"^ — us)' 



where 



u \ u u 



ip^^\u) = -4 ( gi{u) - J dvg2{v) 1 + 2ug2{u) + j dai J da2 

/ 



«3 



+ -2 + V±{ai)) 



«3 



(37) 

a3=l— Qi— a2 



is a combination of twist 4 distribution amplitudes of the pion. The explicit expression for 
(^(4) jg given in Appendix B. 

The twist 4 correction to the light-cone sum rule is easily obtained by taking the 
imaginary part of Eq. (|36| ) in s = {p — q)^ and subtracting the continuum above Sq. After 
Borel transformation one obtains: 



The second term in the r.h.s. of Eq. ( pHD has not been taken into account in It 
appears as a 'surface term' when the correlation function with a denominator [q — upY'^ = 
{—uQ^ + us)"' with > 1 is converted (integrating by parts) into a canonical dispersion 
integral T^u{s,Q'^) = l/vr J dslmT^^{s,Q'^)/{s — s). Adding the expression in Eq. ([38| ) to 
the leading twist 2 contribution ( pT]) one obtains the light-cone sum rule for F.jrlQ'^) to the 
twist 4 accuracy. As seen from Eq. (0), all twist 4 effects have to be identified (to our 
accuracy) as part of the soft contribution to the form factor. Since ^p^^^u) ~ {1 — u) at 
M — >■ 1 (see Appendix B), the twist 4 corrections are of order l/Q"^ in the large limit. 
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Assuming asymptotic expressions for the quark-antiquark-gluon distribution ampli- 



tudes [15, 17, 18] one obtains a compact expression: 



on 

y,(^\u) = =^6\fiyu{3u-2), (39) 

where ^^(IGeV) ~ 0.2 GeV^ is a scale-dependent parameter determining the pion coupling 
to the local quark-antiquark-gluon operator (see Appendix B for the definition). The twist 
4 correction to the sum rule simplifies in this case to: 



revealing at ^ oo the l/Q'^ behavior. Taking in addition the local duality limit — 
oo yields an estimate 

4.2 Factorizable twist 6 contributions 

An estimate of twist 6 contribution to the light-cone sum rule presents a new result of 
this paper. This calculation is interesting for several reasons. As well known [|19[, twist 
4 operators are 'irreducible' in the sense that they cannot be factorized in a product of 
gauge-invariant operators of lower twist. This property is special and limited to twist 4. 
Several light-cone operators of twist 6 exist which can be factorized as a product of two 
gauge-invariant twist 3 operators (or, alternatively, one twist 2 and one twist 4). Sand- 
wiched between vacuum and one-pion state, such operators generally produce two types 
of contributions: Factorizable in terms of a low-twist two-particle distribution amplitude 
times quark (or gluon) condensate, and nonfactorizable that give rise to genuine twist 6 
multiparton pion distribution amplitudes. We emphasize that factorizable contributions 
have to be subtracted in the construction of multiparton distribution amplitudes similar 
as disconnected diagrams proportional to the quark condensate should not be taken into 
account in the nucleon matrix element {N\qq\N) corresponding to the nucleon a-term. 

In the present context, arguments based on conformal symmetry suggest that contri- 
butions of higher Fock states are strongly suppressed at m — >• 1 and their contributions to 
the sum rule are, probably, negligible. Factorizable contributions, on the other hand, are 
expected to supply the missing nonperturbative corrections in the sum rule in the large 

limit and can be large. They are also of principal interest and indicate, as we will see, 
certain limitations for the light-cone sum rule approach. 

Guided by existence of large quark condensate corrections ~ (0|gg|0)2 in classical QCD 
sum rule calculations of the pion form factor [^, in this paper we concentrate on 
factorizable contributions of twist 6 four-quark operators, e.g. 

{0\q{vix)'j^q{v2x)q{v3x)'j^'j5q{v4x)\7i+{p)) = 

= Y^(0|gg|0)(0| ^q{vix)af,^'j5q{v4x) + q{v3x)a^^'j5q{v2x)^ k^(p)) 
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Figure 7: Quark condensate corrections to the correlation function in Eq. (|^). 



IT 



fi'M^(0|gg|0)(0| q{vix)-i^q{vix) - q{y^x)-i^q{v2x) |7r+(p)) 



(42) 



which involve the quark condensate and the two existing two-particle pion distribution 
amplitude of twist 3, (pp and (p^, see Appendix B. Definitions of the both of them include 
the normalization factor 



ml 



(43) 



rriu + rrid 

so that the corresponding contributions to the sum rule appear to be proportional to the 
quark condensate squared. 

We start from the light-cone expansion of the quark propagator (see Appendix A) which 
contains contributions proportional to the covariant derivative D'^Gau of the gluon field 
strength. They are reduced to a quark-antiquark pair due to the QCD equations of motion. 
One quark (antiquark) from this pair can be combined with an antiquark (quark) from the 
initial currents forming a quark condensate as in Fig. 6b, c. 

A straightforward calculation gives: 



vu 



1 1 





^ ^ ^ UV — 3 uv ^ 3q'^v{2 + u) 3q^uv 

3 " \ [q — {l—uv)pf [q — uvpf [q — {l—uv)pf' [q — uvpf 



(44) 



Another source of the factorizable twist 6 contribution is provided by the four-quark opera- 
tors in the light-cone expansion of Eq. (|l|) with a perturbative gluon exchange between two 



currents (see Fig. 7). The technique of this expansion is explained in [16, 22 1. A lengthy 
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but equally straightforward calculation yields: 







u 



1 1 



[q — (1 —uv)pY' 
v{upq — q^) 



[q — {1—uv)pY' 



(45) 



In addition, we have considered the twist 6 parts of the two- and three-particle matrix 
elements corresponding to the diagrams of Fig. 1 and Fig. 5a and have not found any fac- 
torizable contributions. The sum of Eqs. ( PD and (^Sf ) represents, therefore, the complete 
answer for the factorizable twist 6 contributions of four-quark operators. 

The corresponding correction to the light-cone sum rule can be obtained following the 
standard procedure, that is taking imaginary part in s = {p—qY-, subtracting the continuum 
above s = sq in the dispersion integral, and performing the Borel transformation. Due to 
large dimension of the denominators in Eqs. and (|45|), one ends up with a rather 
complicated structure of surface terms at s = Sq. The final answer can be written as 



where 
with 





+h{s,Q^)^^(e{s,-s)e 



-sIM"-" 



(46) 
(47) 



2Q\Q^ + s) 



du 



u- 



1 



1 



2(g2 + s) 



du 



2{2~u)^p{u) + 1 + 



2s 



u- 



2s 



du 



u^ 



1 



U(pp{u) + -(Pa{u) 1 



2uipp{u) - Ir^Paiu) 



US 

uQ"' 



6Q' 



du 



u- 



[2+u)^,{u) 



du 



u- 
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It is easy to see that the expressions in Eqs. (^Sf) and ( ^91) receive contributions from 



both hard u < uq and soft u > uq regions, which we do not write separately in this case. 
The hard contribution takes into account the integration region corresponding to a large 
momentum ~ Q flowing through the gluon line and can be thought of as part of the hard 
mechanism contribution to the form factor generated by product of two twist 3 distribution 
amplitudes, with 'wrong' quark helicities [^, In the light-cone sum rule approach one 



distribution amplitude is present directly, and the second one is modelled using the duality 
approximation, as in the leading twist. 

Inserting the asymptotic expressions for the distribution amplitudes (pp{u) = 1 and 
(Pa{u) = 6u(l — u) and integrating over u one obtains the expansions 

/.(.,g^)^-i, + ^(5-21n^)+..., 

Ms,Q') = --^ + ^Js-ln^)+..., (50) 



g2 g4 

substitution of which in Eq. ( ^6] ) yields the twist 6 correction to the light-cone sum rule to 
the 0{1/Q^) accuracy: 

much smaller than Fi^"*. 

Most importantly, the 0(1/Q^) contributions have cancelled. Inspection of Eqs. ( ^Sf ) 
and ( |49| ) reveals that all 0(1/Q^) contributions in individual diagrams originate from the 
'hard' integration region u < uq and their cancellation involves both diagrams and both 
distribution functions (pa- and ipp, in agreement with 

The observed cancellation of 0(1/ Q'^) corrections is not entirely trivial. One might 
fear that factorization of a local quark condensate brings us back to the deficiency of 
the standard QCD sum rule approach discussed in Sect. 2: Expansion of local operators 



messes up the power counting in the momentum transfer, as observed in [20, 21 1. In other 
language, factorization of the quark condensate is equivalent to using a very bad model 
for the distribution amplitude of the pion created by the interpolation current in Eq. (0), 
corresponding to the sum of two (5-functions, see p6|.p] 

A detailed comparison with the results of Q goes beyond the tasks of this paper. In particular, 
one may ask whether light-cone QCD sum rules can be used to calculate an effective infrared cutoff in the 



hard scattering contribution obtained in [g3| 24 1. To address this issue one has to construct a different 
sum rule, using a chiral-odd interpolating current for the pion. Interpretation of Eq. ( pi] ) in this context 
is difficult because of possible contamination by the Ai meson. We thank M. Beneke and G. Buchalla for 
the discussion which initiated our interest to this problem. 
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Figure 8: The Borel parameter dependence of the hght-cone sum rule for the asymptotic (left) 
and CZ (right) distribution amplitudes, at =1, 3 and 10 GeV^ shown by dashed, solid, and 
dotted curves, respectively. 



5 Numerical analysis 



Combining the twist 2 calculation (l21| ) with the twist 4 corrections in Eq. (|38|) and twist 6 
in Eq. (^), we are now in a position to evaluate the complete light-cone sum rule for the 
pion form factor: 

F,(Q2) =Fi2)(g2,M2) + FW(Q^M2)+Fi6)(Q^M2). (52) 

To avoid misunderstanding, note that terms of higher twist are not suppressed, in general, 
by increasing powers of l/Q^, but rather by increasing powers of the Borel parameter. In 
particular, all twists contribute to l/Q"^ accuracy, with main contributions coming from 
the soft region, in agreement with the general wisdom (see also Sect. 2) that such correc- 
tions come from large transverse distances. The (numerical) hierarchy of contributions of 
different twist is, therefore, a self-consistency check for the light-cone sum rule approach. 

There are several input parameters which should be specified in the sum rule. First 
of all, the pion duality interval, Sq = 0.7 GeV^, is determined by fitting the 2-point sum 
rule (^) to the pion decay constant = 133 MeV. This sum rule is reliable for the 
corresponding Borel parameter M\^^ = 0.7 — 1.0 GeV^. Having in mind that in the light- 
cone sum rule for the same pion channel the Borel parameter should be larger, typically 
of order ~ ^Ipt/ assume 0.8 < < 1.5 GeV^ as a fiducial interval. We have 
checked that changing sq by ±0.1 GeV^ does not produce a significant effect, so that we 



stick to the above standard value 1T3, GO, 211 in what follows. 



The principal input is provided by the leading twist distribution amplitude (see Ap- 
pendix B) 

(pTr{u, fi) = 6uu 1 + a2{fi)C2^'^{u — u) + a4^{^)Cl^^{u — u) + . . . , (53) 

where the coefficients present the main nonperturbative input of interest. Taking into 
account poor accuracy of the present data as well as considerable uncertainties in the 
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Figure 9: The relative importance of contributions of different twist in the hght-cone sum rule 
for the asymptotic (left) and CZ (right) twist 2 distribution amplitudes. Twist 2 (dashed), twist 4 
(dot-dashed) and twist 6 (dotted) contributions and their sum (solid curve) are plotted at = 1 
GeV2. 



sum rules themselves, we cannot aim to distinguish between contributions of different 
Gegenbauer polynomials. We put, therefore, all a„,n = 4,6, .. . to zero and consider the 
values 02 = (asymptotic distribution) and 02(1 GeV) = 2/3 (CZ distribution) as the 
two extreme alternatives. Calculations in this section are done taking into account the 
anomalous dimension of 02 to one- loop accuracy, see Appendix B. Our main goal will be 
to determine 02 from the comparison of the sum rule results with the experimental data. 

The higher-twist distribution amplitudes and relevant parameters represent another set 
of inputs. They are listed in the Appendix B. The uncertainty in higher-twist corrections 
turns out to be sufficiently small and does not influence our final results. 

Finally, one has to specify the renormalization/factorization scale fi. For this numerical 
analysis we use the w-dependent scale /i^ = {1 — u)Q^ + uM"^ for the leading twist 2 
contributions and simply take /i^ = for the soft-region {u > uq) dominated higher- 
twist corrections in the sum rule. The scale dependence is, in fact, rather mild. 

The Borel parameter dependence of the sum rule is shown in Fig. 8 for three different 
values of Q^. As can be seen from this figure, the prediction for the form factor is sufficiently 
stable. 

The relative contributions of different twist to the sum rule are shown as a function of 
in Fig. 9. The twist 4 contribution does not exceed 25-30% of the total result while 

the twist 6 correction is negligibly small. This hierarchy reveals a good convergence of the 

light-cone expansion, at least at > 1 GeV^. For lower values of the higher- twist 

corrections become unstable and the approach breaks down. 

Finally, in Fig. 10 we compare the light-cone sum rule calculation with the available 

experimental data in the interval 1 < < 7 GeV^ taken from p7|, |28|. The dashed and 
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Figure 10: The light-cone sum rule predictions for the pion electromagnetic form factor using 
asymptotic distribution amplitude (dashed), CZ distribution (dotted) and fit to the data (solid). 

dotted curves correspond to the asymptotic and CZ distribution amplitudes, respectively. 
The solid curve presents the best fit, yielding 

4CSR(^ ^ 1 GeV) = 0.12 ± 0.07 1°;°^ (54) 

with = 13.9 for 14 degrees of freedom. The first error comes from the experimental 
uncertainty, whereas the second error corresponds to the variation of the Borel parameter. 
This value for the coefficient 02 translates to the estimate for the characteristic integral: 

/I 
-^^^{u,n^l GeV) = 3.36 ± 0.21 toil ■ (55) 

6 Matching with the NLO perturbative calculation 

As discussed in Sect. 3.3 the light-cone sum rule in the present form does not yet reproduce 
the full perturbative result in the asymptotic limit — > 00. The missing terms correspond 
to higher-order corrections in the light-cone expansion and are difficult to calculate directly. 
Instead, one can make a matching of the light-cone sum rule to the NLO perturbative 
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calculation by following the standard logic of the asymptotic expansion [^]. To this end, 
write the twist 2 contribution to the sum rule defined in Eq. ( plD as a sum of two terms, 
adding and subtracting the leading asymptotic expression at Q"^ ^ oo (see Eq. (|31|), first 
line) 



pertvn; J 27r^ J Q2 J ^ u ' 



Following the arguments of one can prove that higher-order corrections to the sum rule 
must assemble themselves to reproduce the factorized expression 

1 1 

^pert(Q^) ^ jdx jdyip^{x,fx)TH{x,y,Q^,fj.)ip^{y,fj.) (57) 



with 

THix,y,Q^,IJ,'^ 



N,Q''{l-x)il-y) 



l + ^nix,y,Qyf,') + 

TT 



(5^ 



The remainder Fionp(Q^) is suppressed by a power of Q"^ and presents a true nonperturba- 
tive 'higher twist' correction to the usual perturbative result based on collinear factoriza- 
tion. 

Making the substitution (|57D we effectively take into account all higher-order correc- 
tions to the sum rule to 0{1/Q'^) accuracy and neglect such corrections for power sup- 
pressed terms. This procedure tacitly implies that the numerical effect of the replacement 
(0) is more important than of uncalculated (higher-order and higher- twist) corrections to 
-Fn^np(Q^). Such an assumption is natural, but in fact flawed because of potential dou- 
ble counting of perturbative contributions of soft regions. As one signal for this problem, 
one may notice that the perturbative QCD expression suffers from infrared renormalons 
in high orders [^], which have to be cancelled by the corresponding renormalon contri- 
butions to Fionp(Q^)- Using a full resummed expression for -Fport(Q^) together with the 
leading-order expression for Frlonp(Q^) destroys this intricate cancellation and is, therefore, 
not fully consistent theoretically. This is a usual difficulty of making a separation between 
'perturbative' and 'nonperturbative' contributions, which has been discussed in much de- 
tail recently in context of the calculation of power corrections to deep inelastic scattering. 



Drell-Yan processes, event shapes in e+e" annihilation and inclusive B-decays [pO] , pT| . 

An alternative and theoretically better defined possibility is to make a separation be- 
tween soft and hard contributions to the pion form factor with an explicit cutoff, as in 
Section 3, define Fionp(Q^) as the contribution coming from the soft region, and replace 
the 'hard' contribution to the light-cone sum rule by the perturbative expression restricted 
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to the same hard region. A difficulty in this case is that the soft-hard separation in the 
sum rule involves a cutoff in one momentum fraction only and becomes ambiguous when 
applied to the fully factorized expression (0) involving two momentum fractions^. 

In the present paper we consider the ffist possibility because of its relative simplicity. 
We take into account the radiative correction to the hard-scattering kernel []32| , ^ |35[| 



and the complete NLO evolution of the pion distribution amplitude see Appendix B. 
To this accuracy 



where 



and 35 



pert 



-fpert i,Q 



port 



4«(/i2) + ai°(/i 



(59) 



(60) 



2/2 



LO 



NLO 



(/i 



+ 



J TT 



Q 



|a^^(/.^) + (^^) 



l + a^O(/i2) + a^^(/x^) 



NLO 



fc=3 



LO( ,,2\ 



In 



9 

+ 4 



1 + a!f[p^) + aY-'(ti^) \n ^ + 6.58 + 24.99a^{p^) + 21.43(a|-°(fi^))^ 



+ 32.81a^^(/i^) + 32.55(a^^(/i^))^ + 53.37a^^(/i^)a^^(/i 



,LO 



LOi 



.LOi 



(61) 



Note that we do not distinguish between the renormalization and factorization scales. 
The complete expression for the form factor reads, respectively 

F^{Q') = F,UQ') + F^nM^ + ^i'HQ', M') + Fi'\Q\ M') , (62) 

where we have taken into account that twist 4 and twist 6 corrections to the light-cone 
sum rule receive no 1/Q^ contributions to our accuracy. 

For the numerical analysis, we still have to specify the factorization scale. Since, after 
the subtraction of the asymptotic 1/Q^ contribution, the sum rule contribution is domi- 
nated by soft contributions, we choose the fixed scale //^ ~ = 1 GeV^ for simplicity. 
For the perturbative contribution we use 



IGeV' 



(63) 



^2 = ^ ^2^ 

with parameter k in the range 

1/4 < K < 1 . (64) 

Note that with small values of k the scale is almost Q^.j^dependent. Effectively, this choice 
amounts to doing the perturbative expansion to fixed (second) order and not attempting 
a renormalization group resummation. This allows to minimize the problem with double 
counting of infrared regions. 

® A natural solution would be to introduce a cutoff in the transverse quark-antiquark separation rather 
than in the momentum fraction. 
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Figure 11: The nonperturbative correction to the pion form factor (dotted) combined with the 
NLO perturbation theory (sum: solid curve), compared with the pure light-cone sum rule result 
(dashed). Asymptotic pion distribution amplitude is assumed. The grey band shows the scale- 
dependence of the hard-scattering contribution, see text. 

The numerical results are shown in Fig. 11 assuming the asymptotic pion distribution 
amplitude at the scale 1 GeV. The result of the calculation using Eq. (^) and k = 1/2 
is shown by the solid curve with the shaded band corresponding to variation of the scale 
parameter k in the given range. The dotted curve presents the nonperturbative contri- 
bution and the dashed curve is the 'pure' light-cone sum rule calculation with the same 
parameters. The difference between the solid and the dashed curves presents, therefore, 
the net effect of the substitution ([57|). 

The nonperturbative (power suppressed) contribution to the pion form factor shown by 
the dotted curve in Fig. 11 presents considerable interest by itself. It is, obviously, indepen- 
dent on whether the substitution (|57D is used (cf. discussion in the end of Sect. 3.3), and 
turns out to be comfortably small. This smallness may appear to be unexpected after we 
have found large soft (end-point) corrections in Sect. 3, and is due to a strong cancellation 
between the leading order soft contribution to the sum rule (first line in Eq. (p5|)) and the 
large radiative correction (second line in Eq. (^)) corresponding to the sum of soft and 
hard contributions to accuracy. As seen from Eqs. (^), (|30|) the large negative hard 
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contribution ~ l/Q^ plays the most important role in this cancellation. 

Since, according to our analysis, the pion distribution amplitude does not differ sig- 
nificantly from the asymptotic distribution, the theoretical uncertainty in the light-cone 
sum rule calculation of the nonperturbative correction to the pion form factor is domi- 
nated by dependence on the Borel parameter, as illustrated in Fig. 12. With the central 




Figure 12: The light-cone sum rule prediction for the nonperturbative correction to the pion 
form factor. The grey band shows the sensitivity of the result to variation of the Borel parameter 
within O.SGeV^ < < 1.5GeV^. The white line is the calculation for the standard reference 
value K'P = IGeV^ assumed throughout this paper and the dashed curve is the fit (|65|). 

values of parameters, the nonperturbative correction can be parameterized in the region 
1< < 15 GeV^ asQ 

Q'i^nonp(Q') = gV(l-7046 + 1.0662Q2 + 0.0219Q^)2 (65) 

(all numbers in GeV), and the theoretical error (the grey area in Fig. 12) roughly corre- 
sponds to uncertainty in the overall normalization of order ±25%. 

Choosing, as above, a model for the distribution amplitude at the scale 1 GeV as a 
sum of the leading term and the second Gegenbauer polynomial, and fitting the parameter 

^The given parametrisation should not be used for larger values of since it has a wrong asymptotic 
behavior. 
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02(1 GeV) to the data, we find: 



a^ifi = 1 GeV) = -0.06 ± 0.24 ± 0.03 ± 0.03 , 



/d 
— <^Ju,fi = l GeV) = 2.82 ± 0.72 ± 0.09 ± 0.09 . (66) 
u 

The first error comes from the experimental uncertainty, the second error corresponds to 
uncertainty of the nonperturbative contribution (mainly dependence on the Borel parame- 
ter) and the third error is the scale dependence of the NLO perturbative result. Combining 
the two estimates in Eqs. (|5^) and (|66|) and adding the errors in quadrature, we obtain as 
our final result 



a^ifi = IGeV) = 0.1 ±0.1, 



/d 
— !fJu,^i = lGeV) = 3.3±0.3. (67) 
u 

This determination is dominated by the 'pure' light-cone sum rule result in which case 
we included the data points at lower values having higher accuracy. The situation 
will change when sufficiently precise data at > 2 — 3 GeV^ become available. In 
this region the NLO perturbative prediction complemented by the 'higher-twist' power- 
suppressed correction in Eq. ( |65D becomes, from our point of view, a preferable description, 
with potential theoretical accuracy of order 10%. Note that the theoretical status of our 
result for the nonperturbative (soft + hard) correction is similar to model (or sum rule) 
determinations of matrix elements of higher-twist operators in deep inelastic scattering. 



7 Conclusions 

Elaborating on the earlier proposal we have given in this paper a detailed quantitative 
analysis of the pion form factor in the region of intermediate momentum transfers in the 
light-cone sum rule approach and also combining this technique with a complete existing 
NLO perturbative calculation. Our results support the shape of the pion distribution 
amplitude that is close to the asymptotic expression and are inconsistent with the CZ-type 
distributions. Our final estimate for the parameter 02 characterizing the deviation from 
the asymptotic form is given in Eq. ( |67D 

Another important conclusion of our analysis is that the nonperturbative contribution 
to the pion form factor turns out to be rather moderate and does not exceed 30% in 
the full range, see Fig. 11 and Fig. 12. One has to have in mind, however, that 
separation of 'perturbative' and 'nonperturbative' contributions is theoretically not well 



defined because QCD perturbation theory is divergent [§0|- A fully theoretically consistent 
approach necessarily has to introduce an explicit scale separation, and in particular consider 
soft and hard contributions to the pion form factor separately. We have presented a detailed 



^"The smallness of nonasymptotic contributions to (^^ is in agreement with the hght-cone sum rule 
analysis for the 7*77r° transition form factor, compared with the CLEO data |Q. For a recent update 
including NLO effects see ETI. 
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study of the soft-hard separation implemented with a hard momentum fraction cutoff in 
Sect. 3. One finds that soft contributions are generally very large and the smallness of total 
nonperturbative correction is due to cancellations between soft and hard terms of higher 
twist. Thus, somewhat paradoxically, the nonperturbative effects in the pion form factor 
can be small and the soft contributions large, simultaneously! 

To summarize, we believe that the light-cone sum rule approach presents a powerful 
and theoretically consistent framework to the analysis of hard exclusive reactions for in- 
termediate momentum transfers. Main and essential assumption of the method is duality, 
i.e that pion contribution can be isolated from the correlation function by integrating the 
QCD spectral density in the certain energy range - interval of duality. While the numerical 
accuracy of this approximation can be disputed, it satisfies all known QCD constraints and 
provides a perfect laboratory for the study of different interaction mechanisms involving 
several scales. In particular, the scale-dependence of the soft-hard separation studied in 
this work is of general validity. 

Acknowledgments 

The authors thank P. Hoyer and A. Radyushkin for useful discussions. M.M. is grateful 
to the Nordic Institute for Theoretical Physics (NORDITA) for hospitality. A.K. acknowl- 
edges the NORDITA support during his visit when this project was initiated. The work of 
A.K. is supported by the German Federal Ministry for Education and Research (BMBF) 
under contract number 05 7WZ91P (0). 



Here we collect some useful formulae. 

The imaginary part of the radiative correction Hi to the hard scattering amplitude 
reads: 



Appendix A 



-lmHi(Q'^,s,u, iJ,) 
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(Q2 + sfuu 
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where p = Q'^u — us. 



The hght-cone expansion of the quark propagator derived in ([0): 
S{x,Q) = -2(0|r{g(x)g(0)}|0) = 

1 



27r2(-x2)'i/2 i67r2(-x2)<i/2-i 



1 

+2iuu ^XpDxGP'-iux)^ - j du\i{uu - ]^)DpG^''{ux)^, 



+ '-uu{l - 2u)x^ ^D^G^'^iux) + ^uue^uaf^x^D^DxG^^Yl^} + •.• , (A.2) 

where G^"" = gsG^""'{\'' /2), Tr (A'^A^) = 25"^ and d is the space-time dimension. Only 
the terms proportional to the one gluon-field strength and its first covariant derivative are 
shown for brevity, 



Appendix B 

Here we define the light-cone distribution amplitudes of the pion and specify their param- 
eters. The leading twist 2 amplitude (Pn{u) and the twist 4 amplitudes gi{u) and g2{u) 
enter the light-cone expansion of the matrix element 

1 



1 

+fn (^x, - J du e-'^^P-g^iu) . (B.l) 



The QCD equations of motions relate gi and g2 to the quark- antiquark-gluon twist 4 
distributions ip\\, ip±, (^y, and ip±. The latter are defined by the following matrix elements 



{0\d{-x)-f^-f5Gai3{vx)u{x)\n^{p)) = Pf,— — — -U I 'Danp\\{ai)e 



px J 

HgiaPp - 9t^Po)U j P«,y.^(a,)e-^-("'), (B.2) 



29 



+ i9iaPp-9ipPa)U Va.(p^{ai)e-'^^^^^^\ (B.3) 



where Gap = \tai3pxGP^ and the following abbreviations are used: 



r(aj) = tti — + fas, T>ai = da2da2da36 (1 — ai — 02 — as) 



and 



9al3 — 9al3 



XaP(3 + X0P 

a 

px 



The distribution amplitudes are usually constructed using the formalism of the confor- 
mal expansion [^]. To achieve a reasonable accuracy one tries to retain a few first terms 
of this expansion in addition to the leading asymptotic term. The most familiar example 
is the twist 2 pion distribution H 



'^T,{u, jj) = 6uu 1 + a2(/u)C|''^(-u — u) + ai{ij)C'l' {u — u) + . . . 



3/2/ 



:b.4) 



where two orders of the conformal expansion in Gegenbauer polynomials Cn"^ are explicitly 
shown, with 

3 



cT{x) = ''-{hx'-l) 



Cl'\x) = —(21x^- 14x^ + 1). 



15 



(B.5) 



The coefficients a„ determine the nonasymptotic part of Their scale-dependence is 
given in the leading order by 



(B.6) 



where /3o = 11 — and the anomalous dimensions are 



7i°^ = C, 



3 + 



(n + l)(n + 2) 



:b.7) 



In the numerical analysis in this paper we use, in particular, the asymptotic distribution 
(all a„ = 0) and the CZ-distribution (02(1.0 GeV) = 2/3, a„>2 = 0). In next-to-leading 
order the evolution requires an infinite sum of coefficients: 

= 4°(/^')^2(/i') + g2o(/i') 

af 0(/i2) = a^O(/i^)P4(/i') + Q40(/i') + 4°(/i')Q42(/i') 

a^t'^if^') = Q2k o(/x') + a^O(/i2)g2. 2(/i') + a^^(/i')Q2. 4(/i') , A; > 3 (B.8) 
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with the following notations: 

.(1) 



{2k + 3) (n + 2)(n + l) 
{k + l){k + 2) 2(2n + 3) 



7fc 7" 

(0) (0) , n 

7fc - 7n + /5o 



Cli^ = (2n + 3) 



7^°^ -Po + ^CpAkn ^ 2C_F{Akn - Hk + 2) + V^(l)) 



{k - n){k + n + 3) 



(n+ l)(n + 2) 



where 



^;in = h '^^I^^ ) - H^^) + mk -n)- + 2) - V^(l) 



= ^(lnr(^)) f3, = i02-^Np 
= 0, 7^^^ = 111.03, 7^^ = 150.28. 



(B.9) 



(B.IO) 



Expressions for the twist 4 distributions including the next-to-leading corrections in con- 



formal spin have been derived in [p!5| , p!7| : 

5 1 
gi{u) = + -5^e uu{2 + 13^^) 



6 6 

+10u^(2 -3u+ -u^) Inw + 10u^(2 -3u + -u^) Inu 
5 5 



92{u) = —5'^uu{u — u) 



(B.ll) 



V9||(a;i) = 1205^e(Q;i — a2)aia2az , 



- + 2e(l - 203) 



= -1205^Q;ia2a3 



- + e(l-3«3) 



(^_L(ai) = 30(5 03(1 - 03) 



1 



+ 2e(l - 2^3 



(B.12) 
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To this accuracy, the specific combination { ^7\ ) of twist 4 distribution amphtudes reads: 



20 

^^^\u) = —6\^u{3u - 2) - 4:6^eu{2 + llu - 26u^ + 13u^) 



M^(10 - 15u + 6u^) ln(u) + u^{l + 3u + 6u^) ln(l - u) 



(B.13) 



The normahzations of all these distributions are determined by a single nonperturbative 
parameter 6"^ defined as _ 

{7i\gsdGa^.ru\0) = 15'' Uq^ . (B.14) 

The second parameter e in Eqs. ( |B.11| ) and (p.l2| ) is responsible for the first nonasymptotic 
corrections. QCD sum rule estimates yield 5^ ~ 0.2GeV^ B^, [33]: and e ~ 0.5 |1T5| . The 



scale-dependence of these parameters is given by |15 



a 



2\\ 10//3o 



a 



Finally, we should include in our list the twist 3 distribution amplitudes ipp and ip^- used in 
the calculation of the twist 6 corrections. These distributions parameterize the following 
matrix elements: 

1 

{0\d{0)t^5u{x)\n+{p)) = U^i^ jdue-'^'P^^p{u) , 



1 

(0|(i(0)crc,/375M(x)|7r+(j9)) = ^-{PaXp - Pf3Xo)f^l^^ j du 6''''^'' ip ^ (u) , (B.16) 



where /Itt = m1/{mu + rrid). The well known asymptotic form of these distributions: 

ipp{u) = 1, ipa{u) = Quu, (B.17) 

is sufficient for the approximation adopted in this paper. The relation (^31) , together with 
the standard value of the quark condensate (gg)(lGeV) = (— 240MeV)^ yields fi-jrif^ = 
IGeV) ~ 1.56 GeV. Note that the normalization of the twist 6 correction is effectively 
determined by the product ois{[j){qq)''{ij) having in total a negligible anomalous dimension. 
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